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Abstract. We extend Zcilberger's approach to special function identities to 
cases that are not holonomic. The method of creative telescoping is thus 
applied to definite sums or integrals involving Stirling or Bernoulli numbers, 
incomplete Gamma function or polylogarithms, which are not covered by the 
holonomic framework. The basic idea is to take into account the dimension 
of appropriate ideals in Ore algebras. This unifies several earlier extensions 
and provides algorithms for summation and integration in classes that had not 
been accessible to computer algebra before. 



1. Introduction 

In a classical article entitled "A holonomic systems approach to special func- 
tions identities" [16], Doron Zeilberger has shown that the theory of holonomic 
D-modules leads to algorithms for proving identities in large classes of special func- 
tions. In this setting, a function f(x\, . . . , x n ) is represented by a system of linear 
differential equations with polynomial coefficients that annihilate it. The function 
is "holonomic" when it possesses two important properties: (i) besides its defining 
system, / can be specified by a finite number of initial conditions; (ii) the num- 
ber of linearly independent functions among all a;™ 1 • ■ ■ x™ n -^-^ ■ ■ ■ ® 2 n (/) with 

m,\ + • • • + m n + ki + ■ ■ ■ + k n < N grows like 0(N n ). The first property has the 
consequence that many operations on holonomic functions reduce to linear alge- 
bra. It leads to closure properties under sum, product, and specialization. The 
second one is related to the notion of holonomic D-modules and opens the way to 
algorithms for definite integration and summation. For this, Zeilberger developed 
a general method called creative telescoping, for which he gave two algorithms: one 
for the general holonomic case and a faster one in the hypergeometric case. 

Originally, the notion of a holonomic system is only defined for differential sys- 
tems, but there are several ways of extending it to systems of difference (or q- 
difference) equations. Among those, we prefer the use of Ore algebras [6]. There, 
the first property above corresponds to zero-dimensional ideals, which are called 
9-finite. The same closure properties (sum, product, specialization) hold and can 
be performed by Grobner bases computation. Chyzak [5] extended Zeilberger's 
fast hypergeometric creative telescoping to all 9-finite ideals, termination being 
guaranteed inside Zeilberger's holonomic class. 

Another direction of extension concerns functions or sequences that cannot be 
defined by a holonomic system or even a 9-finite ideal. Majewicz [11] has given an 
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algorithm that is able to produce Abel's summation identity 



E ( J *( fc + 0* _1 (n - k + j) n - k = (n + i+j) 



n 



k=0 



automatically and to find similar new identities. Kauers [9] has given a summa- 
tion algorithm applicable to sums involving Stirling numbers and similar sequences 
defined by triangular recurrence equations. This algorithm finds, for instance, 



where 52 and E\ refer to the Stirling numbers of second kind and the Eulerian 
numbers of first kind, respectively. A summation algorithm of Chen and Sun [4] 
is able to discover certain summation identities involving Bernoulli numbers B n or 
similar quantities, for example 



None of the quantities covered by these algorithms admits a definition via a d- 
finite ideal, but all three algorithms are based on principles that resemble those 
employed for holonomic systems and 9-finite ideals. In each case, it turns out that 
the differential/difference equations defining the integrand/summand are of a form 
that permits to prove the existence of at least one non-trivial differential/difference 
equation for the integral/sum by a counting argument. 

In this article, we give algorithms dealing with ideals of Ore algebras that are 
not 9-finitc. They generalize the algorithms known for the 9-finite case and cover 
the extensions to non-holonomic functions discussed so far. Holonomy being lost, it 
is not always the case that creative telescoping can succeed. However, holonomy is 
only a sufficient condition. We show that by considering more generally the dimen- 
sion of the ideals and another quantity that we call polynomial growth, it is possible 
to predict termination of a generalization of Chyzak's generalization of Zcilberger's 
fast algorithm. As special cases, we recover holonomic systems (dimension 0, poly- 
nomial growth 1), but also the special purpose algorithms mentioned above for 
Abel- type sums (dimension 2, polynomial growth 1), Stirling-number identities and 
summation identities about Bernoulli numbers (dimension 1, polynomial growth 1). 
In addition, we get for free a summation/integration algorithm that can deal with 
non-holonomic special functions such as the incomplete Gamma function T(n, z), 
the Hurwitz zeta function £(n, z), polylogarithms Li„(a;), . . . Examples are given 
in Section 4.4. 



Our motivation for using Ore algebras is the convenient polynomial representa- 
tion of linear operators that they offer. Classical notions of commutative polynomial 
rings generalize to this setting. In this section, we recall without proof the basic 
definitions and facts we use (see [6, 8] and their references for proofs, details, and 
history). 

2.1. Ore Algebras. 





2. Ore Algebras and their Ideals 



A NON-HOLONOMIC SYSTEMS APPROACH TO SPECIAL FUNCTION IDENTITIES 3 



Operator d ■ f(x) x ■ f(x) 

Differentiation ? f'( x ) x f( x ) 

Shift S f{x + 1) xf (x) 

Difference A f{x + 1) - f(x) xf(x) 

g-Dilation Q f{l x ) x f( x ) 

Continuous q-difference f(q x ) — ,f( x ) x f( x ) 

^-Differentiation 8^ ^(^fi^ x f( x ) 

g-Shift S (q) f(x + 1) q x f(x) 

Discrete q-difference A' 9 ' f(x + 1) — f(x) q x f(x) 

Eulerian operator xf'(x) x f( x ) 

Mahlerian operator M f( xb ) x f( x ) 

Divided differences /0)-/M xf(x) 

x — a J \ / 



Table 1. Some common Ore operators 



a -derivations. Let A be a commutative algebra over a field fc, and a an injcc- 
tive algebra endomorphism of A that induces the identity on k. A /c-linear en- 
domorphism S of A is called a a -derivation if it satisfies the skew Leibniz rule 
8(uv) = a(u)8(v) + S(u)v for all u and v in A. 

Skew-polynomial rings. The associative ring generated over A by a new indeter- 
minate d and the relations du = a(u)d + 5(u) for all u E A is called a (left) 
skew-polynomial ring. It is denoted by A[d; a, 6]. It does not have zero-divisors. It 
is called an extension of the skew-polynomial ring B[d; er', 5'] when B C A, a' = o\b 
and 8' = 8\b- In examples, we allow ourselves to use different symbols in place of d 
for increased readability. In particular, we use S n ,Sk, etc. for denoting both the 
indeterminate d and the corresponding a in shift algebras. 

Ore operators. A skew polynomial L G A[d; a, S] acts on left A[d; er, (5]-modules — in 
most of our applications, these are modules of functions, power series, or sequences. 
Solving then means finding an element h in such a module such that L ■ h = 0. In 
this perspective, skew polynomials arc called Ore operators. 

For any A € A, the action d : a i— » d ■ a = A<r(a) + 5(a) turns the algebra A 
itself into a left A[d; a, 5]-module. For any u,v in A, one has the product rule 
d ■ uv = a(u)d ■ v + S(u)v. 

Tables 1 and 2 illustrate some common types of Ore operators when A = k[x], 
together with the values of a and 5 that define the associated skew-polynomial ring. 

In all these examples, a and 8 can be written as ^4-linear combinations of d and 
the identity. We call linear a skew-polynomial ring with this property. 

Ore algebras. Ore algebras are a generalization of skew-polynomial rings well suited 
to the manipulation of multivariate special functions. Let C(x) = C[x\ 1 . . . ,x m ) 
be a field of characteristic of rational functions, (ci, . . . , cr„) be n C-algebra mor- 
phisms of C(x) commuting pairwise, and for each i, let 8i be a <7;-derivation, such 
that the S^s commute pairwise and commute with the <jj's when j =/= i. The Ore 
algebra O x = C(x)(d) is the associative C(x)-algebra generated by indeterminates 
d = {di, . . . , d n } modulo the relations 



(1) 



d l a = cri(a)di + S l (a) [a G C(x)), didj = djdi. 
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Operator 


a 


5 




Differentiation 


la 


A 

dx 


IO + 1 


Shift S 


S 





(x + 1)9 


Difference 


s 


A 


(a;+l)9 + l 


g-Dilation 


Q 





qxd 


Cont. g-difference 


Q 


Q-Id 


qxd + (g — l)a; 


g-Differentiation 


Q 


(q-l)x (*3 


- Id) qxd + 1 


g-Shift 


Q 





gx<9 


Discr. g-difference 


Q 


Q - Id 


qxd + (g — l)x 


Eulerian operator 


Id 


dx 


£<9 + x 


Mahlerian operator 


M 





x b d 


Divided differences / i 




X 1— > J » ; 

3: — 


f(a) ad + 1 

a 



Table 2 . Corresponding skew-polynomial rings and their commu- 
tation rules 



An Ore algebra O or an extension A®c(x) © of it is called linear when for each i, 
both &i and can be expressed as C(a;)-lincar combinations of di and the identity. 

Given two Ore algebras O x and ©t, we write O x .t or with an abuse of notation 
C(x,t)(d x ,d t ) for O x ® c Of 

2.2. Ideals. We write / < R to denote that / is a left ideal in the ring R, and 
/ = (pi , . . . , pk) to denote that / is generated by p\ , . . . , pk- The study of a function 
or sequence / translates algebraically into the study of its annihilating ideal ann^ /, 
i.e., the left ideal of operators in an appropriate algebra A that annihilate /. (We 
also write ann / when no ambiguity on A can arise.) Computations concern finding 
generators of this ideal or, at least, of a sufficiently large subideal of it. This is in 
particular the case for creative telescoping that we study here. It computes an ideal 
annihilating the definite integral or sum of interest starting from a description of 
an ideal annihilating the summand or integrand. 

2.3. Grobner Bases. 

Terms. If / = ^2c a d a is a polynomial, each c a d a for which c a ^ is called a 
term of /, c a is its coefficient, d a its monomial, and a its exponent. The total 
degree of / is the maximum \a\ over its terms, where we use the notation 
|a| = | (eti , . . . , a n )\ =«! + •••+ a n . We also use | jS" | for the cardinality of a set S, 
but this should not create confusion. 

Monomial orders. A monomial order is a total order on the monomials that is 
compatible with the product and does not have infinite descending chains. A graded 
order is a monomial order such that d a > d whenever |o;| > |/3|. 

Grobner bases. For a given monomial order, the leading term of a polynomial / 
is the term with largest monomial for that order. We write E(/) for its exponent 
in N n . The crucial property that lets the theory parallel that of the commutative 
case is that E(fg) = E(/)+E(.g). The set Si := {E(/), f e 1} thus has the property 
Si = Si + N" . The complement of Si is a finite union of translates of coordinate 
subspaccs. A Grobner basis G of / is a set of generators / such that Si is the sum 
of E(g) + N" over g G G. 
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2.4. Hilbert Dimension. We write R( s > for the set of polynomials of an Ore 
algebra R of total degree at most s. Then F s > := I D R^ s ' is a vector space over 
C(x). As in the commutative case, the Hilbert function of the ideal / is defined 
by HF/(s) := dim(i?( s ) / 1^) = dimi?( s ^ — dim/( s ); for s large enough this function 
is equal to a polynomial whose degree is called the (Hilbert) dimension of the ideal. 
We denote this integer dim/, or dim^/ when we want to make the ambient ring 
explicit. A reference for the results of this section is [10]. 

Example 1 (d- finite ideals). An ideal is d- finite if its dimension is 0. This special 
class of ideals has been studied a lot from the computational point of view: in this 
case, the quotient R/I is a finite-dimensional vector space, so that many techniques 
of linear algebra apply. 

Example 2 (Hypergeomctric terms). An n-variate sequence u m i,...,m n is a hyper- 
geometric term if 



(mi, . . . ,m„) 



for i = 1, ...,n. If S = (Si,...,S n ), Si representing the shift operator, the 
annihilating ideal of such a sequence in the algebra Q(m)(S) contains opera- 
tors Si — ri(m), . . . , S n —r n (rn). It is therefore 9-finite and moreover, the dimension 
of the quotient as a vector space is 1. 

Example 3 (Stirling numbers). Stirling numbers of the second kind, S^n, fc), 
satisfy the linear recurrence 

£2(71, k) = S2(n — 1, k — 1) + kS2(n — 1, k). 

The ideal generated by this relation in Q(n, k)(S ni Sk) has dimension 1. Properties 
of the generating series imply that this is not a 9-finite sequence, so that 1 is the 
lowest possible dimension to work with. 

The dimension of an ideal can be computed from a Grobner basis for a graded 
order: it is the largest dimension of coordinate subspaces of N™ that belong to the 
complement of the set of exponents of the leading terms [3, p. 449]. 

We note the following inequalities that also hold in this non-commutative con- 
text: 

I C J dim J < dim /, 
/ < C(x, t)(d x ,d t ) => dimc {x>t){am) (I n C(x, t)(d x )) 
< dim c(a , it )( 8a 9t ) /. 

The first one follows from the inclusion of the vector spaces 

lis) c j( s )_ The 

second inequality involves dimensions relative to two different ambient rings. It 
can be seen by considering the following vector spaces: F = C{x 1 i){d x ,dtY'^ 
contains G = C(x,t){d x ) < ^ s ' > and H = Jw. Then the inequality follows from 
G/{G n H) c F/H. 

We also make use of the following properties of dimension. 

Lemma 1. Let I < O x = C(x)(d x ) with dim / = d. 

(1) Foranyd t cd x , \d t \ > d + 1 / H C(x)(d t ) ^ {0}. 

(2) There exists dt C d x of cardinality d such that 

mc(x)(d t ) = {o}. 



(i 
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Proof. We sketch the proof which shows how the notions presented so far interact 
(exactly as in the commutative case). By definition of the Hilbert function, 



The binomial is a polynomial in s of degree \dt\ > d + 1 with positive leading 
coefficient, so that for large enough s, it is larger than HF/(s). The sum of the 
dimensions of the vector subspaces on the left-hand side is therefore larger than the 

(s) 

dimension of Oi and thus they intersect nontrivially. 

The second part follows by considering the exponents of leading terms. Now 
the combinatorial theory is exactly as in [7, Ch. 9]. There exists a coordinate 
subspace of N™ of dimension d in the complement set. This means that there exists 
a subset dt C d x such that no monomial in dt is a leading term of an clement 
of L Thus there cannot be an element of / in those variables only, as was to be 
proved. □ 



Our main result, Thm. 3 in the next section, generalizes the fact that holonomy 
is preserved under definite integration. First, we show how addition, multiplication, 
and action by d behave with respect to dimension, generalizing the corresponding 
closure properties for 9-finite ideals. 

Theorem 1 (Closure Properties). Let <! and let /i,/2 be annihilated by 

I\,l2, respectively. Then: 

(1) dimann((9 ■ /i) < dim/i for all d in {dx} . 

(2) dimann(/i + /a) < max(dim/i, dim/a)- 

(3) If fi, fi belong to the coefficient ring of a linear extension of O x , then 
dim ann(/i fi ) < dim 1\ + dim I2 ■ 

Proof. We show part 3. The arguments for parts 1 and 2 are similar and simpler. 

Setting k := dim I\ + dim J2 + 1 5 it suffices to show that any k elements d\, . . . , dk 
among d are dependent modulo ann(/i/2). 

Given a Grobner basis for I\ with respect to a graded order, each polynomial P G 
O x can be reduced to a normal form P such that P ■ f\ = P ■ f\ and deg P < deg P. 
Moreover, by definition of the dimension, the set of all monomials in {P \ P £ R( s '} 
has cardinality growing like 0(s dlm/l ). The same considerations hold for f 2 . 

By induction on the degree of P, the condition of a linear extension implies 
that P ■ (/1/2) rewrites as a linear combination of monomials (d^ 3 ■ f\){d^ ■ /a), 
with |/3 + 7 1 < deg P. Moreover, we can assume that the monomials have been 
reduced to their normal forms and the inequalities still hold. 

Let s > and consider the following identities 



(|a| < s), where the sums on the right are constructed as above. Then the first 
sum actually ranges over a subset {/3, \f3\ < s} of cardinality 0(s dim/l ) and the 
second one ranges over a subset {7, |— y | < s} of cardinality 0(s dlml2 ). Thus there 



dim/ (s) +dimC*(a;)(dt) (s) 

dimOW + dimC(x)(d t ) {s) - HF 7 (s) 

dimOW + f |9t| +S ) -HFj(s). 




3. Closure Properties 




\P\<sh\<s 
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is a generating set of £)( s dlm7 i s dlm/ 2) _ 0(s k ~ 1 ) monomials for all the summands. 
This implies that for s large enough, there exists a nontrivial linear combination of 
the 0(s k ) polynomials in (2) of the form 

w a (x)d^---d^-(f 1 f 2 ) = 0, 

\a\<s 

as we wanted to show. □ 

Algorithm. The proof gives an algorithm that computes generators of a subideal 
of the desired annihilating ideal, with a dimension that obeys the inequality. For 
increasing s, compute the normal forms of all monomials in d of degree at most s, 
compute linear combinations between them (the kernel of the matrix (u a; ^ i7 )), and 
return these relations if they are sufficiently many to obtain the dimension of the 
theorem. For 9-finite ideals, this returns the same result as the algorithms in [6]. 
Various optimizations are possible. The other closure operations are similar. 

Example 4. The sequence 

fn,m,k,l = (^jS 2 (k,l)S 2 (n - k,m), 

is annihilated by an ideal of dimension at most 2. This follows from Thm. 1 by 
observing: (i) that (?) is hypergeometric, thus 9-finite, and thus annihilated by an 
ideal of dimension 0; (ii) that S 2 (k, I) and S 2 (n — k, m) are Stirling-like (see Section 
4.4) and thus annihilated by certain ideals of dimension 1. More specifically, the 
factors (^J, S 2 (k,l), and S 2 (n — k,m) arc annihilated by the ideals 

({k -n- 1)5„ + n + 1, (k + l)S k +k- n, S m -I, Si- 1), 

(5„-l,5 fc 5, -(Z + 1)5,-1, 5 m -l>, 

(S n Sk - 1, (m + l)S m S k + 5 fc - S m , Si - 1), 

respectively. The algorithm sketched above yields 

I := (1 + n + (1 + m)(l + n)S m - (1 - k + n)S n S m , 

(k - n)S m + (1 + k)S k Si + (1 + fc)(l + m)S k S l S m 

+ (1 + l){k - n)SiS m , 1 + n + (1 + 0(1 + n)S t 

- (1 + k)S k SiS n ) <C(n,m,k,l)(S n ,S m ,Sk,Si) 

as an ideal of annihilators of f n . m , k ,i- It has dimension 2. 

4. Creative Telescoping 

Creative telescoping is basically a combination of differentiation under the inte- 
gral sign and integration by parts, or analogues for other operators. We now give 
it an algebraic interpretation. Our main theorem can be viewed as predicting cases 
when identities arc bound to exist. 

4.1. Telescoping of an Ideal. The heart of the method of creative telescoping 
translates algebraically into the notion of the telescoping of an ideal. 

Definition 1. Let / < O x .t = C(x,t){d x ,dt) be a left ideal. Assume \t\ = \d t \. 
We define the telescoping of I with respect to t = (t\, . . . ,tfc) as the left ideal of 

c{x)(d x ) 

T t (I) := (I + d tl o^ t + ■■■+ d tk o x , t ) n C(x)(d x ). 
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Definition 2. The variables dt = (9t l; ■ ■ ■ ,dt k ) of the Ore algebra <Oa»,t are tele- 
scopable if there exist elements ai, . . . , aft in C(x, t) such that 

5 u {ai) eC{x)\{0} and o u {ai)d tj = d t] a u (ai) (j^i)- 

Note that this is a condition on the algebra, and does not depend on any specific 
ideal. In view of Table 1, this is not a strong restriction for our applications. For 
example, the differential operator d/dt and the difference operator At are telesco- 
pable, with a = t, but the shift operator is not. This notion lets us generalize an 
idea of Wegschaider [14] used at the end of the proof of Thm. 3 below. 

4.2. Polynomial Growth. In Thm. 3 below we give an upper bound for the 
dimension of Tt(I), thus providing a termination criterion for the algorithms in 
Section 5. Our bound depends on the dimension of / and on its "polynomial 
growth", defined as follows. 

Definition 3. The left ideal I < Xi t has polynomial growth p with respect to 
a given graded order if there exists a sequence of polynomials P s {x,t) 1 s G N, 
such that for any a. with \a\ < s, the normal form of P s (x,t)d a with respect to 
a Grobner basis of / for the order has coefficients in C(a;)[i] whose degrees with 
respect to t are 0(s p ). We say that / has polynomial growth p when there is a 
graded order with respect to which it does. 

If all the er's are automorphisms, the polynomial growth is bounded by the 
dimension of the Ore algebra, \x\ + \t\. But the interesting cases are those where 
the polynomial growth is smaller than that. For certain ideals of dimension 0, we 
get a better estimate in Thm. 2 below. For an arbitrary ideal in an arbitrary Ore 
algebra, we do not know how to determine its polynomial growth algorithmically 

yet. 

Example 5. The basis at the end of Example 4 is a Grobner basis with respect to 
a graded order. Inspection of its leading coefficients shows that the sequence 

P s (n,m,k) := JJ (1 + n - k + j)(l + k + jf 
\i\<» 

satisfies the conditions of Definition 3. Since degP s = O(s), it follows that / has 
polynomial growth 1. 

Just before stating our result, we now give notation and a definition for algebras 
amenable to it. 

Definition 4. Set R = C(x)[t] and let R< n be the set of elements A £ R such 
that deg t A < n. A difference- differential algebra is an Ore algebra O x _t such that 
for any i, either 

• Si = and, for any u € R, deg t o~i(u) = deg t u; 

• or o~i = Id and Si is a derivation such that, for any u £ R, deg t Si(u) < 
v + deg t u, 

v G N being fixed. We set S and T> to be the sets of i of, resp., first and second 
types. 

Theorem 2. Let Xi t be a difference- differential algebra endowed with a graded 
ordering. Let I < X) i have dimension 0. Call <fi the map sending any A G O x ,t to 
its normal form modulo I w.r.t. the given graded ordering, and T the (finite) set of 
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P. 



s+1 



monomials in normal form. Then, there exist m,<ES and L £ R<i such that, for 
any i and any (3 £ V , 

(3) mo p ) g i J2 R <^- 

■yer 

Define a sequence (P a )«eN by Pq = I and 

flcm(P. 1 Icm(Lo- i (P 8 )),P,lcm(i,Q a )), tfX>^0, 
\lcm(P s ,lcm(i(r i (P s ))) ) z/P = ; 

where Q s denotes the squarefree part of P s and the 1cm in o~i is over all i. 
If dcg t P s = G(s p ) for some integer p > 0, then I has polynomial growth p. 

Before the proof, note that, by the definition of (f>, the sum in (3) is limited 
to Yf\ < |/3| + 1. Moreover, the existence of a uniform (L,m) is easily obtained in 
the zero-dimensional the set of (i,(3) is then finite. 

Proof. Introduce the sets F s = {d a , \a\ < s} for seN. For n G N, define R< n [T] 
as X^-ygr ^<n9' Y . For a fixed s, suppose there exists P G R and an integer D for 
which 

We look for homologucs P' and .D' for F s +\. As satisfies 

4>(f. +1 ) = 0(p s u|Ja,p s ) c 0(p s )u|J ( /)(a^(p s )), 

we study ^(^A) for A = p- lT Jd p when U G P<d, /3 G T, and |/3| < s. If i G S, 
then from = follows d t A = a l (P)- 1 c7 i (U)d l d' 3 ; therefore, 

<f>(atA) g (za^p)) -1 ^^]. 

Else, i G P and from cr 4 being the identity follows d t A = p- 1 Ud' t d (3 +P- 1 5 l (U)d l3 + 
8i{P- x )Ud p ; thus, <j>{diA) is in 

— i?<D+ m [r] + — R< D+e+l/ d + -p^R<D+D 9 P 

where Q is the square- free part of P and Dq = dcg t Q. Defining P' as the 1cm of 
P, Plcm(L,Q) if V ^ 0, and the a i (P)L's for i e 5 yields A, £ (P') _1 .R[r]. 

Next, setting £>i = deg t (P'/PQ) and Z?2 = dcg t (P'/PP), then A as the maximum 
of max{m, £} + D 2 and max{t + v + D 2 , Do + £>i} if 2? ^ 0, finally £>' = P> + A 
yields 

<f>{F s+1 ) c (P , )- 1 fl< D ,[T]. 

Therefore, introducing the sequence (P s ) defined in the theorem, we get by a first 
induction on s that <fi(F s ) C P^Pp]. Next, under the additional assumption that 
deg t P s = Q(s p ) for some integer p > 0, P and P' in the proof above can be taken 
as P s and P s +i, so that all of £> , -D 1; and D 2 are bounded by degP s+ i — degP s , 
thus by 0(s p ~ 1 ). Therefore, by another induction 

<t>(F s ) cP^RkowFI, 
proving a polynomial growth p for /. □ 
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In the uses of creative telescoping for summation, the identities (or ideals) are 
often stated in terms of shifts, while the operation of interest is a finite difference. 
However, it is possible to compute the polynomial growth in a difference-differential 
algebra with shifts, and it will be the same when considered in an algebra with 
difference operators. This is the meaning of the following theorem, whose proof 
based on transporting graded ordcrings we omit. 

Proposition 1. For a given difference- differential algebra O x ,t and indeterminates 
u = (iti, . . . , Uk), consider 

(V(t,-u) = @ x ,t ®c C(u)[d Ul ;S Ul ,0] . . . [d Uk ;S Uk ,0] and 

<°'x,(t, u ) = <0«,t ®c C(u)[d Ul ; S U1 , A U J . . . [d Uk ; S Uk , A„J, 

where A Ui = S Ui — Id (i = 1, . . . , k). Let /i be the (bijective) left-C(x,t,u) -linear 
map sending each d Ui to d Ui + 1 and all the d x , dt to themselves. Then an ideal 
^x,(t,u) an d mOO f3 ^ x (t u ) have the same polynomial growth. 

4.3. Main Result. Our main result is the following sufficient condition for creative 
telescoping. For simplicity, we state it for \t\ = \dt\- 

Theorem 3 (Creative Telescoping). Let I < Q x .t = C(x 1 t)(d Xl dt) be an ideal of 
dimension d and polynomial growth p. If \t\ = \dt\ and the dt's are telescopable 
then 

dimo x T t {I) < d+ (p— 1)1*1, 
whenever this bound is nonnegative. In particular, when the bound is smaller than 
\x\, the ideal Tt(I) is non-trivial. 

Proof. If d + (p — l)\t\ > \x\, there is nothing to show. Otherwise, by Lemma 1, it 
is sufficient to show that any k := d + (p — l)\t \ + 1 elements di, . . . , dk in {d x } are 
dependent modulo Tt(I). 

Let s > and consider the following set of members of I from the definition of 
the polynomial growth of /: 

(4) p^^-w-E^*)^ ' a i^ s - 

l/3|<s 

The coefficients c a .p(x, i) can be viewed as C(a;)-linear combinations of monomials 
in t of degree 0(s p ). The ideal / having dimension d, these sums in (4) for all a 
such that \a\ < s actually range over a common subset of {/3, \f3\ < s} of cardinality 
0(s d ). Thus, there is a generating set of 0(s d+p '*') monomials in t, d for all the 
summands. This implies that for s large enough, there exists a nontrivial linear 
combination of the 0(s fc+ '*') polynomials in (4) of the form 

P s (x,t) J2 C a (x)d? ■ --d^d t at =: P(s)Q e /. 

Multiplying by 1/P s (x, t) shows that Q € In C(x)(d x , dt). Now the operators 
commute with the coefficients of Q. Thus, Q can be rewritten 

(5) Q = R + d tl Q 1 + --- + d t]tl Q ltl , 

with R e C(x)(di, ...,d k ) G T t {I). HR^Q, we have found the element of T t {I) 
we were looking for. Otherwise, since Q ^ 0, there exists a Qi ^ 0. Since <9 ti is 
telescopable, there is an element o, G C(x,t) such that 

o~ti {a.i)d tl = hi + d u ai, a u (a l )d tj = d t] a ti (a,) 
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for some bi £ C(x) \ {0}. Multiplying (5) by cr t . (<Zj) yields 

I 3 0u(ai)Q = hQi + d ti Qi H h d tw Q\ t \- 

Now, 7^ can be rewritten as in (5). Repeating this process if necessary, we 
eventually reach a non-zero element of Tt(I), as was to be proved. □ 

4.4. Examples. 

Proper hypergeometric sequences. These arc sequences with two indices that are 
hypergeometric as in Example 2 with the further constraint that they can be written 

(6) ^ = Q(^)|^^^ 

lli=i ( u i m + v l k + w l y. 

where £ € C, Q is a polynomial and the a^'s, 6, 's, u,:'s, i>i's are integers. A typical 
example is the binomial coefficient (™) . Since such a sequence is hypergeometric, 
a Grobner basis of annii m j; for any order is formed by the relations 

S m — u„ i+ i ! fe/u mi fc and Sk — u m ^+i/u m ,k- 

The normal form of S^J S^ 2 with respect to this basis is simply the rational func- 
tion u m+sl: k+ S2 1 u m.k- Since the dj, 6j, Uj, Uj's are all integers, the common denom- 
inator of all these rational functions for si + s 2 < s has a degree that grows only 
linearly with s [12], as it is bounded by 

P s (m, k) = Q(m, k) ~\{uim + u^fc + w t + j) 

\j\<(\ Ui \+\vi\)s * 
x ^[(dim + bik + Ci + j). 

\j\<(\ai\ + \bi\)s * 

In our terminology, the corresponding annihilating ideal has polynomial growth 1, 
by Thm. 2 and Prop. 1. Further generalizations to the multivariate and g-cases can 
also be considered. Thm. 3 (with d = 0,p= 1, |i| = l,n = 2) generalizes the result 
that creative telescoping applies to proper hypergeometric sequences. 

General hypergeometric sequences. Not every hypergeometric sequence is proper. 
For example, 

1 (1m -2k- V 



mk + 1 \ m — 1 

is not. By a criterion of Abramov [1], creative telescoping fails on this example, 
i.e., Tt(I) = {0}. This phenomenon is well consistent with our theorem, because 
the nonlinear factor in the denominator implies that the annihilating ideal I has 
polynomial growth 2: a possible choice for P s is 

P s (m,k)= [ ((m + i)(k + j) + 1) x ' [ (m-2k + i), 

i+j<s \i\<3s 

whose degree is quadratic in s, hence the polynomial growth by Thm. 2 and Prop. 1. 
Thm. 3 (with d = 0, p = 2, \t\ = l,n = 2) implies the trivial bound dim Tt(/) < 1, 
which is reached in this example. 
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Holonomic Functions. This notion was popularized for special functions in [16]. 
The technical definition is related to holonomic D-modulcs. Thm. 2 implies that it 
is sufficient to consider 9-finite ideals, as was shown first by Takayama [13]. 

Proposition 2. Let I < O x .t = C(x,t)(d x ^} be a left ideal of dimension of the 
differential Ore algebra (o~i = Id, Si = d/dxi). Then I has polynomial growth 1. 

Proof. The crucial property is that in the differential setting, the leading coefficient 
of an clement does not change upon multiplication by some 9, (see Table 2). The 
result is a corollary of Thm. 2, which is applicable since I has dimension 0. The 
polynomials P s defined in the theorem here take the form P s = L s , having degree 
0(s 1 ) and thereby proving a polynomial growth 1. □ 

Put together, Prop. 2 and Thm. 3 imply dimTt(7) = 0, and thus recover the 
celebrated closure of holonomic functions under definite integration as a special 
case. 

Stirling-like sequences. This notion was introduced in [9]. A sequence is called 
Stirling-like if it is annihilated by an ideal of operators 

(u + vSfrS^ + wS2Sf llll s 2 + t 2 S m2 ,. . . ,s n + t„S m J, 

where u, v, w, s 2 , ■ ■ ■ , s n , t 2l ■ ■ ■ , t n are polynomials that split into integer-linear fac- 
tors, and a,/3,7, S are integers subject to a certain nondegeneracy condition. A 
typical example is (™)S2(fc, m), where S 2 refers to the Stirling numbers of the sec- 
ond kind. It was shown in [9] that such ideals lead to non-trivial creative telescoping 
relations whenever n > 2. These ideals have dimension 1 and polynomial growth 1, 
so our theorem (with d = l,p = l,\t\ = l,n > 3) includes this result as a special 
case. 

Abel-type sequences. This notion was introduced in [11]. A sequence is called Abel- 
type if it can be written in the form 

u m ,k{k + r) k {m -k + s) m " k -^-— 

k + r 

for some proper hypergeometric term u mj fc. The annihilating ideals of such se- 
quences can be written 

{aS m Sk - bS r , cS m - dS s ) 

for certain polynomials a, 6, c, d. It was shown in [11] that such ideals lead to non- 
trivial creative telescoping relations. These ideals have dimension 2, and if u m .k is 
as in (6), then 

P s (m, r, s, k) = P s (m, k) x (k + r + j)(m — k + j), 

b'l<« 

with P s (m, k) being the polynomial sequence stated above for proper hypergeomet- 
ric terms, justifies that their polynomial growth is 1, so our theorem (with d = 2, 
p = 1, \t\ = 1, n = 4) includes this result as a special case, too. 
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Bernoulli examples. Chen and Sun [4] do not give a formal description of the class 
of summands to which their algorithm is applicable. Most of their examples concern 
sums with summands of the form hk,rn 1 .m 2 Bak+bm 1 +cm2 where h is a hypergeometric 
term, B refers to Bernoulli numbers, Bernoulli polynomials, Eulcr numbers, or Euler 
polynomials, and a, b, c are specific integers. The annihilating ideals of such objects 
have dimension 1, and, if h is proper hypergeometric, polynomial growth 1. In this 
case, our theorem guarantees the success of creative telescoping. 

Further examples. Our theorem also extends to sequences and functions for which 
no special purpose summation or integration algorithm has been formulated so far. 
For instance, the sequence f n ,m,k,i considered in Examples 4 and 5 is not Stirling- 
like, but only annihilated by an ideal of dimension 2 and polynomial growth 1. Our 
Theorem predicts the existence of creative telescoping relations (cf. also Examples 6 
and 7 below). We list some further identities from the literature that were previ- 
ously considered inaccessible to computer algebra, but that can be proven with 
creative telescoping. In all the following identities, the integrand is not holonomic, 
but annihilated by an ideal of dimension 1 and polynomial growth 1, so our theorem 
predicts a priori that relations for the integral must exist. 



In these identities, B refers to the Beta-function, ( to the Hurwitz zeta function, 
Li„ to the nth polylogarithm, T to the incomplete Gamma function, and J n to the 
nth Bessel function. Greek letters refer to parameters, n, m, k are discrete variables, 
x,y, z are continuous ones. Note that the polylogarithm, while holonomic for each 
specific n, is not even 9-finite when n is "symbolic" . Also, the integrand of the last 
integral, despite being elementary, is not holonomic. 

All these identities are proven by computing operators annihilating both sides, 
by making use of closure properties for sums, products or d (Section 3) and by 
creative telescoping in the case of definite sums and integrals (Section 5 below). If 
the operators found in this way generate an ideal of dimension 0, then proving the 
identity reduces to verifying a finite number of initial values. For higher dimensional 
ideals, the number of initial values to be checked may be infinite. For example, for 
completing the proof of the first integral identity above, it remains to check the 
identity for k = 1 and all n > 0: 



an identity not much easier than the original one. In many instances, however, the 
identities to be verified as initial conditions are trivial. 




fc)C(n — k, a), 
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5. Algorithms for Summation and Integration 

5.1. A Fasenmyer-Style Algorithm. The first algorithmic approach to sym- 
bolic summation goes back to Fasenmyer and was formulated for the summation 
of hypergeometric terms [12]. The following generalized version of her algorithm is 
applicable to any given annihilating ideal in an Ore algebra. 

The algorithm follows from the proof of Thm. 3. Given an ideal I<C (x, t) (d x ,dt) } 
the algorithm searches for elements of / n C(x)(d x ,dt) (called i-free operators). 
Proceeding by increasing total degree, it makes an ansatz with undetermined coef- 
ficients for such an operator, reduces it to normal form with respect to a Grobncr 
basis for /, brings the normal form to a common denominator, and then compares 
the coefficients with respect to the £'s of the numerator to zero. This gives a lin- 
ear system of equations over C{x) for the undetermined coefficients, which is then 
solved. If the system has no solution, the procedure is repeated with a larger total 
degree. If a solution is found, it leads to an operator that can be written in the 
form (5) 

A + d tl B 1 + -<- + d tw B\q el 

where ^ A e T t (I), by following the steps at the end of the proof that ensure A ^ 

0. The procedure can be repeated to search for further operators in T t (I) until all 
operators found generate an ideal of the dimension predicted by Thm. 3. 

Example 6. For the ideal / from Example 4, with k as summation variable, this 
algorithm discovers a first operator for total degree 4 that may be written in the 
form 

S m + Si + (2 + I + m)SiS m - SiS m S n 

+ (S k - l)((m + l)S m Si - S m S n Si + Si). 

As a consequence, we find 

A := S m +S t + (2 + l + m)SiS m - SiS m S n e T k (I). 
Therefore A belongs to the annihilating ideal of the sum 

^ (j^jS 2 (k, l)S 2 {n - k,m). 

This proves that the sum is equal to ( l+ l m )S2(n,l + m): A also belongs to the 
annihilating ideal of that quantity, and it agrees with the sum for n = and 
arbitrary I and m. 

5.2. A Zeilberger-Style Algorithm. Zeilbcrgcr's "fast algorithm" was originally 
formulated for summation of proper hypergeometric terms only [15], then a differen- 
tial analog was given for hyperexponential functions by Almkvist and Zeilberger [2] . 
This was later extended by Chyzak [5] to an integration/summation algorithm for 
arbitrary 9-finite functions, with termination guaranteed in the holonomic case. 
We extend this algorithm further to the case of integrands or summands denned 
by arbitrary annihilating ideals. 

The fast approach is applicable only for single sums or integrals, i.e., if |t| = 

1. (Refer however to [5, Sec. 3.3] for an iterated treatment of summations and 
integrations.) Let I<\C{x, t)(d x ,dt) be given by a Grobner basis G with respect to a 
graded order. The algorithm searches for operators A+dtB G / with A e C(x)(d x ) 
and B £ C(x, t)(d x ,dt). Proceeding by increasing total degree, it makes an ansatz 
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for both A and B and computes the normal form of Q := A + dtB with respect 
to G. Without loss of generality, only irreducible terms need to be included in the 
ansatz for B. In order for Q to belong to /, it is necessary and sufficient that all 
the coefficients in the normal form of Q be zero. Comparing them to zero leads to a 
system of first order linear functional equations (a "coupled system" ) , which is then 
solved. While in the <9-finite case this system is always square, in the case of positive 
dimension it may be rectangular, owing to extraneous equations potentially being 
introduced by irreducible terms in dtB of degree d+1. However, it is always possible 
to separate the equations into a square coupled system (that can be solved in a 
first step as in Chyzak's algorithm) and a system with additional linear algebraic 
constraints (that can be accommodated in a second step by computing suitable 
linear combinations of the solutions obtained in the first step) . Any solution A + 
dtB E I found in this way gives rise to an element A 6 T t (I). The procedure 
may be repeated until the elements of T t (I) found in this way generate an ideal in 
C(x)(d x ) whose dimension matches the dimension predicted by Thm. 3. 

Example 7. Applying the algorithm to the ideal / of Example 4 with respect to k, 
we obtain a coupled system of size 5 x 14 when A and B arc assumed to have total 
degree 2 and 1, respectively. This system has no solution. For total degrees 3 and 2, 
respectively, the obtained system is of size 14 x 28 and has the nontrivial solution 

A = S m + Si + (2 + I + m)SiS m - SiS m S n , 
B = fc(fc+l) + (m + l)fe 

fc 2 — 1 — n — kn k — n — 1 ' 
hence A + (S k - l)B E I, and A G T k (I). 

6. Final Comments 

We have seen that algorithms for special functions doing computations in Ore 
algebras are not restricted to special functions with 9-finite annihilating ideals. 
Instead, both closure properties algorithms and algorithms for definite integration 
and summation can be formulated for arbitrary ideals. The treatment could be 
extended further by including, for instance, ideals of Laurent Ore algebras. 

Our generalized algorithms rely on the notion of ideal dimension as well as on 
the notion of polynomial growth we introduced in Definition 3. According to this 
definition, the polynomial growth depends on the monomial order imposed on the 
underlying algebra. Future research will focus on reducing the notion of polynomial 
growth to an intrinsic property of the ideal at hand, as well as to devising an 
algorithm for computing the polynomial growth. 
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